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Abstract-Texture analysis and classification are some of the 
problems which have been paid much attention by image 
processing scientists since late 80s. If texture classification is done 
correctly and accurately, it can be used in many cases such as 
object tracking, visual pattern recognition, and face recognition. 
Since now, so many methods are offered to solve this problem. 
Nearly all these methods have tried to extract and define features 
to separate different labels of textures accurately. This article has 
proposed an approach which has an overall process on the 
images of textures based on Local binary pattern, Gray Level 
Co-occurrence matrix, and edge detection to determine energy 
variations. Finally, extracting the features of the images would 
classify them. Although this approach is a general one, it could be 
used in different applications. In the final part a proposed 
approach has been tested on the stone texture and the results 
have been compared with some of the previous approaches to 
prove the quality of proposed approach. Low time and 
computation complexity, noise insensitivity, and rotation 
invariant are some of the advantages of the proposed approach. 

Keywords-Texture Classification; Gray Level Co-occurrence; 
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I. INTRODUCTION 

Since 1990s, with the progress of computer vision and 
image processing, texture classification turned out to be one of 
the main subjects in the literature of these sciences. As texture 
classification has a close relationship with sciences such as 
machine learning, it functions in areas like pattern recognition, 
object tracking, defect detection, face recognition, etc. The 
main problem in texture classification relates to two topics: 
finding and introducing the best features for texture 
description; and selecting the best and the most corresponding 
kinds of classifiers with selected features. 

Regarding to these two topics, we have witnessed different 
kinds of discussions. The techniques used on texture analysis 
or classification, are discussed in four categories, namely 
statistical approaches, structural approaches, filter based 
methods, and model based approaches. Table I shows a 
summary list of some of the key texture analysis methods that 
have been applied to texture classification or segmentation. 
Clearly, statistical and filter based approaches have been very 
popular. 

Statistical texture analysis methods measure the spatial 
distribution of pixel values [3]. They are well rooted in the 
computer vision world and have been extensively applied to 
various tasks. A large number of statistical texture features 
have been proposed, ranging from first order statistics to 
higher order statistics. Amongst many, histogram statistics 
[16], co-occurrence matrices, autocorrelation, and local binary 
patterns have been applied to texture analysis or classification. 
In structural approaches, texture is characterized by texture 



primitives or texture elements, and the spatial arrangement of 
these primitives. Thus the primary goals of structural 
approaches are firstly to extract texture primitives [2], and 
secondly to model or generalize the spatial placement rules. 
The texture primitive can be as simple as individual pixels, a 
region with uniform gray levels, or line segments. The 
placement rules can be obtained through modeling geometric 
relationships between primitives or learning statistical 
properties from texture primitives. The filter based techniques 
largely share a common characteristic, which is applying filter 
banks on the image and computing the energy of the filter 
responses. The methods can be divided into spatial domain, 
frequency domain, and joint spatial/spatial-frequency domain 
techniques. Model based methods include, among many others, 
fractal models, autoregressive models, random field models, 
the epitome model, and the texem model. 

TABLE I INEXHAUSTIVE LIST OFTEXTURAL ANALYSIS 
METHODS 





Method 


Statistical 


1. Histogram Properties 

2. Go -occurrence Matrix 

3. Local Binary Pattern 

4. Other Gray Level Statistics 

5. Autocorrelation 

6. Registration-based 


Structural 


1. Primitive Measurement 

2. Edge Features 

3. Skeleton Representation 

4. Morphological Operations 


Filter Based 


1. Spatial Domain Filtering 

2. Frequency Domain Analysis 

3. Joint Spatial/Spatial-Frequency 


Model Based 


1. Fractal Models 

2. Random Field Model 

3. Texem Model 



Energy amount is one of statistical measures that have 
most meaningful information of texture. The energy amount 
variation of a texture image in different situations may provide 
good features for classification. So, in this paper an approach 
is proposed to texture classification based on energy variation. 
The proposed approach includes two stages. In the first stage, 
the texture is processed by LBP, GLCM and edge operators. 
After that the energy amount of output images is computed. In 
the second stage, the energy amounts of first stage are 
compared by energy amount of original image to computation 
energy variation. So in this stage, an accurate feature vector is 
provided to texture classification. In the result section, to 
prove this claim, by gathering some images of textures in four 
models of stones, and applying proposed approach on them, 
stone texture classification has been done. The high accuracy 
in texture classification in the results shows the quality of 
offered approach. 
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The reminder of this paper is organized as follows: Section 
II is related to the description of local binary pattern (LBP) 
algorithm and the way of this estimation. Section III is related 
to the description of Gray level Co-occurrence matrix (GLCM) 
algorithm. Section IV has a brief description about the 
different aspects of the edge detection filters. Section V 
describes the main method of this paper and finally, the results 
and conclusion are given. 

II. LOCAL BINARY PATTERN 

One of the most popular texture analysis operators is Local 
Binary pattern (LBP) that was introduced in [7]. It is a gray- 
scale invariant texture measure computed from the analysis of 
a 3x3 local neighborhood over a central pixel. The LBP is 
based on a binary code describing the local texture pattern. 
This code is built by the thresholding of a local neighborhood 
of the gray value of its center. 

The eight neighbors are labeled using a binary code {0,1} 
obtained by comparing their values to the central pixel value. 
If the tested gray value is below the gray value of the central 
pixel, then it is labeled 0, otherwise it is assigned the value 1: 






ifKX,Vi) < I(x ,y ) 
Otherwise 



(1) 



di' is the obtained binary code, di is the original pixel 
value at position i and dO is the central pixel value. With this 
technique there is 256 (28) possible patterns (or texture units). 

The obtained value is then multiplied by weights given to 
the corresponding pixels. The weight is given by the value 2i- 
1. Summing the obtained values gives the measure of the LBP: 



j _ V8 j' ? i 

l LBP — ^i = l u i^ 



(2) 



Fig. 1 shows an example on how to compute LBP. The 
original 3x3 neighborhood is given in Fig. 1 (a). The central 
pixel value is used as a threshold in order to assign a binary 
value to its neighbors. Fig. 1 (b) shows the result of 
thresholding of the 3x3 neighborhood. The obtained values are 
multiplied by their corresponding weights. The weights kernel 
is given in Fig. 1 (c). The result is given in Fig. 1 (d). 
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LBP = 1+8+128 = 137 



Fig. 1 Computation LBP 



m 



The sum of the resulting values gives the LBP measure 
(137). The central pixel is replaced by the obtained value. A 
new LBP image is constructed by processing each pixel and 
its 3x3 neighbors in the original image. Fig. 2 shows an 
example of the resulting LBP images. 





(b) 
Fig. 2 Processing input image: (a) original image of travertine stone; 
(b) LBP image 

III. GRAY-LEVEL CO-OCCURRENCE MATRIX 

Gray Level Co-occurrence Matrix (GLCM) was proposed 
in [6] by Haralick and Shanmugam. It is very useful in texture 
analysis. It calculates the second order statistics related to 
image properties by considering the spatial relationship of 
pixels. GLCM depicts how often different combinations of 
gray levels co-occur in an image. The GLCM is created by 
calculating how often a pixel with the intensity value i occurs 
in a specific spatial relationship to a pixel with the value j. The 
spatial Relationship can be specified in different ways, the 
default one is between a pixel and its immediate neighbor to 
its right. However we can specify this relationship with 
different offsets and angles. The pixel at position (i,j) in 
GLCM is the sum of the number of times the (i ,j) relationship 
occurs in the image. 

Fig. 3 describes how to compute the GLCM. It shows an 
image and its corresponding co-occurrence matrix using the 
default pixels spatial relationship (offset = +1 in x direction). 
For the pair (2,1) (pixel 2 followed at its right by pixel 1), it is 
found 2 times in the image, then the GLCM image will have 2 
as a value in the position corresponding to Ii = 1 and Ij =2. The 
GLCM matrix is a 256x256 matrix; Ii and Ij are the intensity 
values for an 8b it image. The GLCM can be computed for the 
eight directions around the pixel of interest (Fig. 4). Summing 
results from different directions lead to the isotropic GLCM 
and help achieve a rotation invariant GLCM (Fig. 5). More 
details are in [4]. 
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Fig. 3 Description of the gray level co-occurrence matrix 
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Fig. 4 Directions used for computing isotropic GLCM 
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Fig. 5 Processing GLCM: (a) original image of travertine stone; 
(b) GLCM image processed in 0' 

IV. EDGE DETECTION 

Edge detection is the process of localizing pixel intensity 
transitions. The edge detection has been used by object 
recognition, target tracking, segmentation, and etc. Therefore, 
the edge detection is one of the most important parts of image 
processing. There mainly exist several edge detection methods 
(Sobel [11], Prewitt [12], Roberts [13]). These methods have 
been proposed for detecting transitions in images. Early 
methods determined the best gradient operator to detect sharp 
intensity variations [14]. Commonly used method for 
detecting edges is applying derivative operators on images. 
Derivative based approaches can be categorized into two 
groups, namely the first and the second order derivative 
methods. The first order derivative based techniques depend 
on computing the gradient several directions and combining 
the result of each gradient. The value of the gradient 
magnitude and orientation is estimated using two 
differentiation masks [15]. 

In this work, Sobel, an edge detection method is 
considered. This method is preferred and will be used in this 
work. The Sobel edge detector uses two masks: one vertical 
and one horizontal. These masks generally use 3x3 matrices. 
The matrices specifically having 3x3 dimensions are used in 
MATLAB (see edge.m). 

V. PROPOSED APPROACH 

The approach, which has been proposed in this article, is 
in a main boxconsisting of three sub-boxes as shown in Fig. 6. 

In this approach, first of all, according to the algorithm 
described in Section II , the LBP image would be estimated 
for the input image. The energy amount of the LBP's image is 
then calculated by using Equation (3). 
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Fig. 6 The flowchart of proposed approach 
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(3) 



Where, "n" and "m" are the size of image and N (i, j) is the 
normalized processed image. 

Next, according to the algorithm mentioned in Section III , 
the second order statistics which is called GLCM would be 
estimated for the original texture image. GLCM can be 
estimated in eight directions and achieve one image in every 
direction. Afterward, the energy amount of the GLCM images 
would be computed by using Equation (3). In the last part of 
processing, Sobel filter, the edges of entering original texture 
image would be recognized. In the same way like the previous 
steps, the energy amount of the edge detected image is 
computed. In this point, every image has a feature vector like 
F. It is shown Equation (4). 



F=<Energy (LBP), Energy (GLCM), Energy (Edge)> 



(4) 



In the second stage, the energy amount of original texture 
is computed. Next, it is compared by every dimensions of 
feature vector F. So, a feature vector likevdould provide 
which has meaningful information on input image can be used 
for classification of this kind of textures, as it is shown in 
Equation (5). 

F'=< |O.E-P.E (LBP)|, p.E-P.E (GLCM)|, |O.E-P.E (Edge)| > 

(5) 

Where, O.E. means the original image's energy and P.E. 
means the energy amount of processed image by operators 
which is described in parenthesis. 

VI. RESULTS 

As it was mentioned, the main aim of this article is to offer 
a new approach for texture classification. Thus to get results 
and observe the efficiency of proposed approach, three kinds 
of stones named granite, travertine and hatchet, 60 images 
were taken by digital camera. It means 20 images of every 
model of stones. The images have been processed by the 
approach offered in this article, and the proposed features 
where computed for them according to Section V . To reduce 
the complexity and dataset's dimensions only GLCM was 
estimated in the 0' and 180'. Then there has been a dataset 
made consisting of 60 instance and 4 attributes. Every 
instance has a label that is the model of that stone. Finally, by 
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using some of the classifiers such as KNN, NaiveBayes, and 
J48-Tree, and by using N.Fold method, the accuracy of the 
dataset that has been made for stone textures is computed, 
which is shown in Table II . In the second and the third rows 
in the table, the stone textures were classified based on LBP 
and GLCM images to compare the proposed approach with 
some of the previous approaches in term of accuracy. As it is 
shown in Table I , the accuracy of the new approach is much 
higher than other previous approaches. The main advantages 
of the proposed approach in this article can be mentioned in 
two points: corresponding with nearly all of the classifiers; 
and introducing a series of new features which are able to be 
computed for different applications. 

TABLE II ACCURACY RATE OF CLASSIFICATION DATA SET 



NClassifier 
Approach ^v 


3NN 


Naive Bayes 


J48- 
Tree 


Energy Variation 


94.2 ±0.6 


93 ±0.1 


91.1 ±0.4 


LBP 


93.6 ±0.1 


92.1 ±0.4 


88.2±0.2 


GLCM 


87.3 ±0.5 


82.2 ±0.6 


78.4±0.5 



VII. CONCLUSION AND FUTURE WORK 

The main aim of this paper was to propose an approach for 
texture classification. In this respect, the Sections II , III and 
IV have algorithms of LBP, GLCM and edge estimations for 
the images and in Section V the proposed approach has been 
offered and described thoroughly. As it was described in this 
approach, first the input images were processed based on the 
main box, shown in Fig. 6, then, the energy amounts were 
computed for the exiting images. In continuation, the energy 
variation is computed and the basic feature vector is provided. 
Next, in the result part is collated a dataset for three kinds of 
stones. Finally by using some classifiers and the N,Fold 
method, classification is done. So the results showed that the 
features which have computed by proposed approach can use 
for classifying all kinds of stones by great accuracy. 

Results show that the proposed approach has a low 
computation and time complexity. Also, according to the 
computation global image energy, the proposed feature vector 
is not sensitive to noise. On the other hand, the proposed 
approach is rotation invariant because the energy amount of 
texture is independent of rotation. 

One interesting future research direction is to use the 
proposed approach to texture segmentation or object tracking 
by analyzing the object or image via feature vector F total. 
The future research direction is also to use other operators 
such as Wavelet filters. 
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